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S. K. Samaddar\ J. N. De\ X. Vinas^, and M. CentelW 

^Saha Institute of Nuclear Physics, 

1/AF Bidhannagar, Kolkata 700064, India 

"^Departament d'Estructura i Constituents de la Materia 

and Institut de Ciencies del Cosmos, Facultat de Fisica, 

Universitat de Barcelona, Diagonal 647, 08028 Barcelona, Spain 

^ : Abstract 

(<— ^ I The density and excitation energy dependence of symmetry energy and symmetry free energy 

CN I for finite nuclei are calculated microscopically in a micro canonical framework taking into account 

Qh' thermal and expansion effects. A finite-range momentum and density dependent two-body effective 

interaction is employed for this purpose. The role of mass, isospin and equation of state (EoS) 
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on these quantities is also investigated; our calculated results are in consonance with the available 
experimental data. 
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I. INTRODUCTION 

The symmetry energy is a measure of the energy involved in converting the excess neu- 
trons to protons in asymmetric nuclear matter. A kinetic contribution to it comes from the 
associated shift of the neutron (n) and proton (p) Fermi energies, another contribution comes 
from the difference between the (n-p) interaction and that between like pairs (n-n or p-p). 
Traditionally, the symmetry energy per nucleon or the symmetry energy coefficient Ce of 
infinite nuclear matter has been determined from fits of experimental binding energies with 
various versions of the liquid drop formula [l| . But it refers only at the saturation density 
and at a temperature T = 0. Its value is usually taken to be between 30 and 35 MeV. 

Understanding the details of the structure, mass and the cooling of neutron stars [2| 
or simulating the dynamics of supernovae collapse [3| entails a knowledge of the density 
and temperature dependence of the symmetry energy. The abundance of relatively heavier 
elements in explosive nucleosynthesis or even the existence of exotic neutron or proton-rich 
nuclei produced in collisions of radioactive nuclei have a direct lineage to this knowledge. 
The neutron skin thickness of heavier nuclei has also been found to be intimately correlated 
to the density derivative of the symmetry energy [^, |5|, la] as it reflects the pressure difference 
on the neutrons and protons. 

Collisions between nuclei at relativistic energies offer the best hope of studying properties 
related to isospin asymmetry (symmetry energy, symmetry free energy, etc.) of nuclear mat- 
ter at supranormal densities. Inference can be made there from comparison of theoretical 
prediction with experimental data on symmetry energy-sensitive observables like differential 
flow of neutrons and protons or from the multiplicity ratio of 7r^/7r+, K^ / K^ ^ etc. [7|, |8|, but 
no flrm conclusions can yet be made since the experimental isospin-sensitive signals cannot 
be considered very deflnitive ^]. At subnormal densities, studies on nuclear multifragmen- 
tation offer a unique tool to determine the characteristics of the nuclear symmetry energy 
or symmetry free energy as a function of density and excitation energy. In intermediate 
energy heavy ion collisions, a hot dilute nuclear system is formed which expands to reach 
the equilibrium state and ultimately fragments into many pieces. The produced fragments 
bear signatures of the properties of the hot expanded system prior to fragmentation. These 
include the excitation energy dependence of temperature (caloric curve) and density as well 
as the symmetry energy and symmetry free energy at subnormal densities produced at dif- 



ferent excitations. Data related to isotopic distributions [10], isospin diffusion |lll. Il2l. Il3 



and isoscaling [Ij, Il5| have recently been analyzed and estimates of symmetry coefficients 
at different densities and excitations have been obtained. These estimates give somewhat 
different predictions and are also not fully conclusive. 

There have been numerous studies on the symmetry energy of nuclear matter based on 
the different many-body theories using various nucleon-nucleon interactions or interaction 
Lagrangians [7|]. These studies provide very useful tools for understanding the properties of 
hot and dense nuclear matter. It has been noticed that the calculated density dependence of 
the symmetry energy coefficient differs appreciably depending on the choice of the theoretical 
models and interactions. For the symmetry free energy of infinite nuclear matter, there are 



recent investigations done in the mean field framework [16j. In Refs. |17l. |18|. the symmetry 
energy and symmetry entropy of very dilute nuclear matter have been calculated exploiting 
virial expansion techniques where clusterization of light fragments is taken into account. 

For finite systems, however, there are fewer available calculations for the symmetry en- 
ergy or the symmetry free energy and for their dependence on density and energy. In Ref. 



[19], symmetry free energy coefficients of fragments produced in nuclear multifragmentation 
have been calculated from the variance of the isotopic distributions obtained in a statistical 
multifragmentation model. The present authors have performed a calculation 20|] of the 
symmetry energy coefficient of finite nuclei based on the finite temperature Thomas- Fermi 
(FTTF) formulation. This calculation was done microscopically in a microcanonical frame- 
work using a finite range, momentum, and density dependent effective interaction [2l|. The 
calculated symmetry energy coefficients at different excitations and densities were compared 
with the available scant experimental data. There is an ongoing discussion regarding whether 
the experimental data for the symmetry coefficients should be connected to the symmetry 
energy or to the symmetry free energy 22]. We take the viewpoint that they refer to the 
symmetry free energy, in accord with recent studies 18|, |l9j In the present work we calculate 
the symmetry free energy coefficient for a number of nuclei in the FTTF formulation. The 
calculation of the symmetry energy coefficient in Ref. 20| was done in the local density 
approximation (LDA). Calculations with some improvement over the LDA are reported in 
the present paper. In addition, the dependence of the symmetry coefficients on the mass 
and isospin content of the nucleus as well as on the underlying EoS are considered. 

The organization of the paper is as follows. In Sec. 11, outlines of the model used in 
the calculation are presented. Section 111 contains the results and discussions. Concluding 
remarks are given in Sec. IV. 



II. THEORETICAL FRAMEWORK 

The methodology employed to calculate the symmetry energy and symmetry free energy 
coefficients as a function of excitation energy or density is outlined in the following. 

A. Modeling the hot nucleus 



When two nuclei collide at intermediate energy a hot nuclear system of neutrons and pro- 
tons is formed, which is assumed to be in thermodynamic equilibrium and can be described 
by a temperature T. The density profile of this hot system is generated in the finite temper- 
ature Thomas-Fermi (FTTF) approximation with a chosen two-body effective interaction. 
The details of the employed FTTF procedure are already documented in Ref. |23| and we 
do not present them here. 

For an expanding system pursuing the equilibrium co nfig uration (as described later). 



the surface diffuseness is likely to play an important role [2J]; thus, a zero-range force like 



the Skyrme interaction widely used to explore nuclear ground-state properties may not be 
very suitable for generating such a density profile. It is further noted that a constrained 
expanded system in the FTTF approach may lead to numerical instabilities 25|, |26| and 
the gradient (surface) terms in the energy density functional were replaced with a suitable 
Yukawa interaction. We have therefore chosen a modified Seyler-Blanchard (SBM) effective 
interaction for the FTTF calculations. This interaction is of finite range with momentum 
and density dependence and is given by |2l(] 



Ves{ri,T2,P,p) = -Cl^u 



l-^-rfMp(ri)+p(r2)r 



exp(— r/a) 
(r/a) 



(1) 



An effective isospin dependence in the interaction is brought through the different strength 
parameters Ci for hke-pair {n — n, p — p) and Cu for unhke pair {n — p). The relative 
separations of the nucleons in configuration and momentum space are given by r = |ri — r2| 
and p = |pi — P2|- The densities at the sites of the two interacting nucleons are given 
by p(ri) and p(r2). The parameter a corresponds to the range of the interaction, b and 
d determine its momentum and density dependence; the density exponent n controls the 
stiffness of the nuclear EoS. This interaction reproduces quite well the ground-state binding 
energies, root-mean square charge radii and isoscalar giant monopole resonance energies 
for a host of even-even nuclei. With a density exponent n = 1/6, the incompressibility 
of symmetric nuclear matter Koo is 238 MeV. A stiff EoS with K^o = 380 MeV can be 
simulated with n = 4/3. 

In the FTTF approach, the nucleon density profile at temperature T has the form 

Pr{r) = AUr)Ji/2{vAr)), (2) 

where 

AUr) = ^[2mrAr)Tf\ (3) 

and Jk{i]t) is the Fermi integral 

JK{Vr)= — 7 zdx, (4) 

Jo l + exp{x-r]r) 
with the fugacity rj^ given as 

Vrir) = [fir-VAr)]/T. (5) 

In Eqs. (2)-(5), r is the isospin index, rrir^k the effective fc— mass of the nucleon coming 
from the momentum dependence of the interaction, /i^ the chemical potentials and Vrir) 
the effective single-particle (SP) potential (Coulomb included). 

When r^ -C 0, the system is very dilute with V ~ and then p ~ e^^'^ , a constant. At 
large distances, the particle density therefore does not vanish. The pressure at the surface is 
then nonzero making the system thermodynamically unstable; the density then depends on 
the size of the box in which the FTTF calculations are performed. This problem is overcome 
in the subtraction procedure 27|, l28| , where the hot nucleus, assumed to be a thermalized 



system in equilibrium with a surrounding gas representing evaporated nucleons, is separated 
from the embedding environment. The method is based on the existence of two solutions 
to the FTTF equations, one corresponding to the liquid phase with the surrounding gas 
[Ig) and other corresponding to the gas {g) phase. The density profile of the hot nucleus in 
thermodynamic equilibrium is given by pr = Prig — pT,g- It is independent of the box size 
in which calculations are done. It also goes to zero at large distances, implying a vanishing 
surface pressure. We call this the base density (it is also sometimes called the liquid profile). 
The conservation of the nucleon number of each species A^,- of the hot nucleus gives 

J[Pr,igir)-p,,,ir)]dr = N,. (6) 

The energy E of the required nucleus is given by 

E = Eig- Eg, (7) 



where Eig and Eg are the total energies of the hquid-gas system and of the gas alone. 
The total entropy in the Landau quasi-particle approximation is 

S = -Y1 [9Aer,T) [frlnfr + (1 - fr) ln(l - fr)] dSr, (8) 

where /,- is the single-particle occupancy function 

frier, fir,T) = [1 + exp{ (^^ - /i.)/T}]-\ (9) 

and Qr is the subtracted single-particle level density. Once the energy and entropy are 
known, the free energy is calculated from F = E — TS. 

In the above, the description of the hot nucleus is grand canonical obtained from the 
minimization of the grand potential (the temperature is a constant). In experimental con- 
ditions, however, when two nuclei collide, the hot system is formed in isolation, its total 
excitation energy remains a constant. The system might be compressed initially, resulting 
in a collective flow in the decompression stage but we ignore it in the present work. The 
system is microcanonical; to attain equilibrium, it expands in quest of maximum entropy. 
It is, however, still possible to describe the system statistically by an effective temperature 
T. It has the operational advantage that it helps in defining an occupation function that 
can be employed in evaluating various observables like energy, entropy, etc. 

The expansion of the hot nucleus is simulated through a self-similar scaling approximation 
for the density, 

px{r) = AV(Ar), (10) 

where the scaling parameter A is unity for the unbloated nucleus and decreases with expan- 
sion, lying in the range < A < 1; pa('") is the scaled density and p{r) is the base density 
profile generated in the subtracted FTTF framework. Besides its simplicity, there is no a 
j»non justification for this choice, however, it has been shown that with a harmonic oscillator 
potential, at relatively small temperatures, the scaled density profiles and those generated 
self-consistently in a constrained Thomas-Fermi 25| procedure are equivalent 



One further needs to account properly for the effect of collectivity, as the coupling of 
the single-particle motion with the collective degrees of freedom [30l] is not included in the 
FTTF procedure. This coupling introduces an extra energy dependence in the nucleon 
effective mass {rrii^, the cj— mass) in addition to the fc— mass. The total effective mass m* 
can then be written as 

m =m . (Ill 

m m 



The cj— mass is surface-peaked and has values generally larger [31[ than the nucleon mass 
m. This increase brings down the excited states from higher energy to lower energy near 
the Fermi surface, thus increasing the many-body density of states at low excitations. The 
system can then accommodate comparatively more entropy at a given excitation energy. The 
coupling of collectivity with the nucleonic single-particle motion may thus have a significant 
role in getting the equilibrium maximal entropy configuration. A self-consistent evaluation 
of rriuj is very involved; for simplicity, we take the same phenomenological form of Refs. 



|32l . I33l . I34j | for it. An in-depth presentation of our computational method of the expanded 
hot nucleus with inclusion of collectivity can be found in Ref. |29| and thus we do not dwell 
further on it here. 



B. Symmetry energy 

The symmetry energy esym of nuclear matter characterizes how the energy rises as one 
moves away from equal numbers of neutrons and protons. For asymmetric nuclear matter 
at density p = Pn + Pp with asymmetry parameter X = {pn — Pp)/p, the symmetry energy 
is defined as 

esymip, T, X) = e(p, T, X) - e(p, T, X = 0), (12) 

where e is the total energy per nucleon of nuclear matter, given as 



e(p) 



h: 



2 



2m 



:T{p) + Slip) 



i, (13) 



In the above equation, the first and second terms within the square brackets are the kinetic 
and potential energy densities for infinite nuclear matter at a density p. 
The symmetry energy can be written as 

esym{p,T,X) = Ce{p,T)X^ + 0{X^). (14) 

The terms beyond X"^ are negligible for values of X one encounters in nuclei. The nuclear 



matter symmetry energy coefficient Ce is obtained from [17 

1 & 
2dX 

The symmetry free energy coefficient Cp can similarly be defined as 



Ce{p, T) = -T^^esymip, T, X) \x=o. (15) 



CFip,T) = l^/,^^(p,T,X)U=o, (16) 

where fsymip, T, X) is the symmetry free energy per nucleon defined in the same manner as 
in Eq.(12) with e replaced by /. 

To compute the coefficients Ce and Cp in finite nuclei we adopt the following prescription. 
Once the neutron and proton equilibrium density profiles of a nucleus with A^o neutrons and 
Zq protons {Aq = Nq + Zq) at an excitation energy E* and temperature T are known, the 
symmetry energy coefficient can be calculated in the local density approximation (LDA) 
as t23 

Ce{E*) (^^^^)' = \J p{r)CE{p{r),T)X\T)dr. (17) 

Here, CE{p{'r),T) is the symmetry energy coefficient at temperature T of infinite nu- 
clear matter at a density equal to the local density p{r) of the nucleus and X{r) = 
{pnir) — Ppir)) I pij") is the local isospin asymmetry. One can obtain analogously the sym- 
metry free energy coefficient Cf{E*) of a finite nucleus. 

In the LDA, the particles at each point in space feel the potential as if it were locally a 
constant. The neutron and proton potentials in the configuration space are calculated at a 
temperature T for infinite matter at a value of the local density p{r) to evaluate CE{p{r),T) 
or Cf{p{t),T). In a finite nucleus, these potentials at any point should also contain infor- 
mation on the densities at nearby points, which in the extended Thomas-Fermi (ETF) [35| 



method is taken into account by recasting the kinetic energy density as a functional of not 
only the local density but also its derivatives. The correction to the energy density at a 



temperature T, up to second order in h is J35|, l36 



^2(p)=-^2(p)+T(T2(p), 



(18) 



where J^2ip) and o"2(p) are the corrections to the free energy density and entropy density, 
respectively. They are given as 
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(19) 



(20) 



In the above two equations, p refers to the local density p{r) and / = m/m*{r) is a functional 
of p. The quantity C(^)) to a good approximation, is 



av) 



36 L 



1 + 2/vaT^ 



and 



z/(r/) 



, Ji/2{ri) dC 
J~i/2{r]) dr]' 



(21) 



(22) 



The corrections ^2(p) and ^2(p) are added to the local energy and free energy densities 
perturbatively, in the spirit of variational Wigner-Kirkwood theory |37|, to calculate the 
improved symmetry energy and symmetry free energy coefficients. 



C. Isotopic scaling and symmetry free energy 



It has been observed by various experimental groups 38|, |39|, |40|, |4l|, |42] that the logarithm 
of the ratio R defined as 



R 



Y2{N,Z)/Yi{N,Z), 



(23) 



where Yi and I2 are the yields of a particular fragment with A^ neutrons and Z protons 
from two different fragmenting sources differing in the neutron-proton ratio at the same 
temperature follow a relation of the type 



In R = K + {aN + I3Z). 



(24) 



This observation is known as isoscaling; the coefficients a and (3 are the parameters charac- 
terizing the isoscaling behavior and K is the normalization factor. 



The parameter a has been related to the symmetry coefficient C through the relation 



where the suffixes 1 and 2 correspond to the two fragmenting systems. The quantities 
(Zq, A^i denote the values in the fragmenting system from whose disassembly the fragment 
f Z, A) is produced. Various authors have derived Eq. (l25l) under different approximations 
(id . I12I . |42] and the coefficient C has commonly been related to the symmetry energy co- 



efficient Ce- In this interpretation, the isospin dependence of entropy has been neglected, 
which may be a fair approximation at around normal density but may not be so for low 
densities as encountered in the tail region of the density proffie of a nucleus at a relatively 
high temperature. In some recent literature [l8|, [l9|], the need to include the asymmetry 



dependence of entropy has been stressed and the symmetry coefficient in Eq. fl25p has been 
interpreted as that pertaining to the symmetry free energy. Furthermore, whether the sym- 
metry coefficient refers to the fragmenting source or to the primary fragments at freeze-out 



is not fully settled. In Ref. [42|, it is interpreted as the symmetry coefficient of the pri- 
mary fragments. In Ref. [lO|, the basic interpretation is the same, but the properties of the 
fragments are conjectured to be modified due to "in-medium" effects because of presence of 
other neighbouring fragments in the freeze-out volume. In the sequential Weisskopf model 
in the grand canonical limit [12] as applied for an expanding emitting nucleus, the symmetry 
coefficient is linked to that of the fragmenting source. In our present communication, we 
take the symmetry coefficient to be the symmetry free energy of the expanded mononuclear 
system in its most probable configuration at a fixed excitation energy i5*. 



III. RESULTS AND DISCUSSIONS 

A. Infinite nuclear matter 

The SBM interaction, as has been mentioned earlier, reproduces well the bulk properties 
of nuclei. For symmetric nuclear matter as well as for neutron matter, the EoS obtained [43|, 



44l | with this interaction also compares very favorably with those calculated microscopically 
with realistic interactions in a variational approach [45|, |46|. In Fig. 1, we display the 
symmetry energy coefficient of nuclear matter (at T = 0) as a function of density. Since our 
focus of interest is in the subnuclear density region, the results are presented up to p ~ 0.75po 
where po is the saturation density taken as 0.154 fm~^, its value for the SBM interaction. The 
calculated results are seen to be well within the range obtained in microscopic calculations 
[4 7( 1 with different bare (Argonne vl8) and effective (SLy4 and Gogny) interactions. 

Up to the saturation density the symmetry energy coefficient calculated with the SBM 
interaction can be very well represented by 

C^(p)^C^(po)(-^y, (26) 

with Ce{po) = 34.0 MeV and 7 = 0.65. Though the experimentally extracted value of 
the exponent 7 is still fraught with some uncertainties, significant constraints on it have 
been determined from different observables in recent years. Comparison of results from the 
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transport model with recent experimental data on isospin diffusion constrain the value of 
7 to around 0.69 — 1.05 at subnuclear densities [7[]. The neutron and proton transverse 
emission ratio measurements [48| present some new constraints on 7 somewhat larger than 
0.5, whereas measurements from isotopic distributions 15| provide a value of 7 close to 



0.69. Consideration of the giant dipole resonance properties of ^osp]^ puts a constraint 
23.3 < Ce{p ~ O.lfm-^) < 24.9 MeV (49J, which implies a value of 7 ~ 0.55. 

In Fig. 2, the symmetry coefficients Ce and Cp as a function of density of nuclear matter 
are shown at T = 10 MeV in the upper panel. The difference between Ce [Eq. (15)] and 
Cp [Eq. (16)] is amplified with decrease in density, in consonance with that obtained in 
Ref. |16|]. This is understandable from entropy considerations. Our calculations have been 
done in the mean-field model, inclusion of cluster formation at low densities would increase 



the values of these coefficients somewhat [17[. We find that the equilibrium density {i.e. 



the state at zero pressure) of nuclear matter falls off linearly with excitation energy and 
may be very well represented by p = po(l ~ Q.QAE* / A)., with E* expressed in MeV. The 
symmetry coefficients Cp and Cp for different excitations at equilibrium densities calculated 
using Eqs. (15) and (16) are shown in the lower panel of the figure. The dependence of these 
coefficients with excitation is found to be nearly linear and may be well represented as 
Cp{E*/A) ~ Ci5(0)(l - 0.024^*/^) and Cp{E*/A) ~ Cp{Q){l - 0.028^*/^). We have 
studied the role on Cp and Cp of using a soft EoS (n=l/6, i^oo=238 MeV) and a hard EoS 
(ri=4/3, -ft'oo =380 MeV). The effect of the EoS on both the coefficients in infinite nuclear 
matter is found to be small except at very low densities. 

B. Finite nuclei 

We have calculated the symmetry coefficients for a number of nuclei in order to study 
their mass and asymmetry (Xq = (A^o " Zq)/Aq) dependence as a function of density and 
excitation energy. For the mass dependence, we have chosen ^^^Au and "'"S, both having 
practically the same Xq. For the isospin dependence, we have considered the isobar pair 
isogj^ and ^^°Cs. The relevant experimental data on the symmetry coefficients are very few; 



they are available mostly in the mass region Aq ~ 100 — 120 [15|, [50]. We have therefore 
studied the nucleus ^^°Sn to have a comparison of the calculated results with the experimental 
data. The model is tested further in a wider perspective; we calculate the evolution with 
excitation of temperature (caloric curve) and density of this nucleus as experimental data 



5ll . |52| are available around this mass number. 



1. Grand canonical approach 

In Fig. 3, the caloric curve, the central density pc in units of the ground-state central 
density pc,o and the symmetry coefficients Cp and Cp for the nucleus ^^°Sn are displayed as 
a function of excitation energy E*/A in panels (a), (b), and (c), respectively. These calcu- 
lations have been performed with the base density profile generated in the grand canonical 
framework where all the excitation energy has been locked in the thermal mode, i.e., there is 
no expansion energy. The experimental data for the caloric curve and densities correspond 
to medium-heavy nuclei (100 < Aq < 140). They have been taken from Ref. [511 for the 
caloric curve and from Ref. [52] for the densities. We have also included in the figure the 
available experimental data for the symmetry coefficients; the open triangles and the filled 



circles are from Ref. 15| and the open and filled squares are from Ref. 50|- As discussed 
earlier, we interpret these data as pertaining to the symmetry free energy coefficient. The 



data from Ref. [15| correspond to collisions between mass-symmetric nuclei with total mass 



Aq = 116. The source size was taken there to be somewhat less, Aq ~ 100, because of the 



reduction due to preequilibrium emission. The data in Ref. [50| were extracted for collisions 
of ^^C on ^^^'^^^Sn. The symmetry coefficients there are given as a function of temperature. 
We have expressed them as a function of excitation energy using the Fermi-gas expression 
E* = aT"^ with an effective level density parameter a = A/ 10. 

The present calculations are done for the fragmenting source ^^°Sn to give an orientation 
on the excitation energy dependence of the symmetry coefficients. As reported later, the 
symmetry coefficients are found to be weakly dependent on the mass of the fragmenting 
system, but somewhat sensitive to its Nq/Zq ratio. The dotted black line and dot-dash 
blue line of Fig. 3 correspond to calculations in the LDA for Ce and Cp, respectively. The 
corresponding calculations with second order corrections incorporated are represented by the 
magenta dash line and the full red line. The calculated caloric curve matches very well with 
the experimental data except at high excitations. The calculated densities are, however, 
overestimated. Correlation of the symmetry coefficients with the density is displayed in 
Fig. 4. The notations used for the different lines are the same as in the bottom panel of 
Fig. 3. In panel (c) of Fig. 3, it is seen that the calculated symmetry free energy coefficients 
follow the experimental trend rather well, but in Fig. 4 the mismatch between theory and 
experiment [l5| becomes very apparent indicating the limitations of the grand canonical 
approach. 

2. Microcanonical approach 

The fact that the hot nuclear system formed in energetic nuclear collisions is an isolated 
system and the limitations of the grand canonical approach exposed in Fig. 4, motivate one to 
describe the evolution of the system in microcanonical thermodynamics. In this framework, 
the system expands in search of the maximal entropy configuration. In panels (a) and (b) of 
Fig. 5, the caloric curve and the evolution of the density with excitation energy calculated 
in the microcanonical approach are displayed for the system ^^°Sn. The comparison of the 
observables with the experimental data is now improved, showing the importance of the 
proper treatment of the expansion phase for the equilibrium configuration. The bottom 
panel displays the symmetry coefficients Ce and Cp in the LDA and also with the inclusion 
of second-order corrections. The different lines have the same meaning as in the bottom 
panel of Fig. 3. With increase in excitation, the importance of the second-order corrections 
is found to decrease; this is attributed to the slower fall of the density for nuclei bloated 
with excitation. At higher excitations, the system becomes more expanded and dilute and 
a possible enhancement of the symmetry coefficients with respect to the present calculation 



may come from clustering at the surface [17|, [18 



The correlation of the symmetry coefficients with density is displayed in Fig. 6 for the 
same system ^^°Sn. The notations used for the calculated results are the same as in Fig. 4. 
Allowing for the uncertainties in the experimental extraction of the density and of the 
symmetry coefficients, it is found that the calculated correlation follows the experimental 
trend well. A noticeable improvement of the results over those depicted in Fig. 4 is observed. 

The dependence of the symmetry coefficients Ce and Cp for finite nuclei on the EoS of 
the underlying nuclear interaction is displayed in Fig. 7 at different excitations. We have 
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chosen ^^°Sn as the representative system. All the calculations presented in this figure and 
in Fig. 8 are done with the inclusion of the second order corrections. At the same excitation, 
both Ce and Cp are larger for the stiffer EoS. This is understood from the fact that at the 



same excitation, the equilibrium configuration is more compact for the stiffer EoS [29|. As 
a whole, the symmetry coefficients are found to be not too sensitive to the choice of the EoS 
we have made. 

The excitation energy dependence of the symmetry free energy coefficient Cp for all 
the five nuclei studied is displayed in panel (a) of Fig. 8. The lines from top to bottom 
correspond to the systems ^^''Sn, ^^^Au, ^^°Sm, ^^°Cs and ^''S, respectively. The mass and 
asymmetry dependence of the symmetry coefficient can be easily inferred from the figure. 
The comparison of the results for the systems ^^^Au and ^°S (having practically the same 
asymmetry) indicates the lowering of the symmetry coefficients with decreasing mass. The 
lighter nucleus has a lesser value of Cp because of the predominance of the surface effects. 
Similarly, the isospin or asymmetry dependence can be inferred from the comparison of 
results of the isobar pair ^^°Sm and ^^°Cs. The lower values of the symmetry coefficient 
Cp for the more asymmetric nucleus ^^°Cs can be traced down to the fact that isobars 
with higher asymmetry have effectively softer EoS 29|. It is seen that the results for the 
symmetry coefficients for the pair ^^^Au and ^^''Sm are practically indistinguishable. This 
refiects an interplay of the effects due to mass and asymmetry. This is further amplified in 
the larger values of Cp for ^^°Sn, which has an appreciably smaller mass than ^^^Au but 
has also a very small asymmetry Xq = 0.09. In panel (b) of Fig. 8, the symmetry free 
energy coefficients of all the nuclei studied are displayed as function of their equilibrium 
densities corresponding to different excitations. The variations of the density correlation of 
the symmetry coefficients with mass and isospin for the nuclei studied are very similar to 
those seen for the excitation energy in the upper panel of the figure. The results for the 
symmetry energy coefficient Cp exhibit nearly the same trends with excitation and density 
as Cp in the present figure, and therefore we do not display them. 

As seen in Fig. 8, the excitation energy dependence of Cp for all the five nuclei discussed 
is almost linear and the results corresponding to each nucleus run nearly parallel. As in 
the case of nuclear matter, this dependence can be well approximated by a linear relation 
Cp{E*/A) = Cp{0){l-apE*/A) with ap ~ 0.054 MeV'^ The same holds for the symmetry 
energy coefficient Cp (not shown in the figure), for which we find Cp{E*/A) = Cp{0){l — 
apE* /A) with ap ~ 0.064 MeV~^. The faster fall-off of the symmetry coefficients of finite 
nuclei with increasing excitation compared to those of nuclear matter is attributed to the 
comparatively lower equilibrium density of the isolated nuclei at the same excitation. 

The density dependence of the symmetry free energy coefficient of the finite nuclei can 
be fitted with a general expression of the form 

where k^ and Kg are the volume and surface constants contributing to the symmetry coef- 
ficient. The exponents 71 and 72 depict the density dependence of the volume and surface 
contributions, respectively. We have included a term KsymXl, with Xq being the asymmetry 
parameter of the finite nucleus. This is done in order to test the size of an eventual departure 
of the symmetry energy in finite systems from the quadratic dependence on the asymmetry 
parameter which is assumed in the definition of the symmetry coefficients. 

A least-squares fit of Eq. fl27j) to the calculated values for the symmetry free energy 
coefficient fixing the values of n^ (34 MeV) and 71 (0.65) to those of infinite nuclear matter, 
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and considering all the five systems studied in the excitation energy range 1 < E*/A < 10 
MeV, gives Ks = 1-46, 72 = 0.17 and Ksym = 1-55 with a root mean square deviation ^ 6%. 
The significantly lower value of 72 compared to the volume exponent 71 points to a weaker 
surface density dependence. A free variation of all the five parameters improves the least- 
squares fit very little compared to the variation of three parameters mentioned above. In 
infinite nuclear matter, the symmetry energy and symmetry free energy are known to be 
well represented with a quadratic term in the asymmetry parameter X, the quartic term 
being negligible. In finite nuclei, the existence of surface and Coulomb effects may modify 
this scenario. The result Ksym = 1-55 found indicates, however, that the effect is relatively 
small for the typical values of Xq in nuclei. 



IV. CONCLUDING REMARKS 

We have investigated the energy and density dependence of the symmetry energy and 
symmetry free energy coefficients of finite and infinite nuclear systems. The dependence of 
these coefficients on the EoS, mass and isospin content of nuclei have further been explored. 
The calculations are done in a microscopic microcanonical framework using a momentum and 
density dependent finite range effective interaction. The density dependence of the symmetry 
energy coefficient of infinite nuclear matter calculated with this interaction compares very 
well with those obtained from other microscopic calculations. 

Our main focus in the present work is to explore the density and energy dependence of 
the symmetry coefficients of finite nuclei. First, we have investigated the predictions of our 
considered model for these coefficients for infinite nuclear matter in the subnuclear density 
range. In the density range 0.1 < p/po < 1, the symmetry energy coefficient for nuclear mat- 
ter at T = is found to be well reproduced by Ce{p) — Ce{po) {p/PoY with Ce{po) = 34.0 
MeV at the saturation density and 7 ~ 0.65, well within the experimental range of values. 
For finite nuclei, the calculations have been performed in the local density approximation 
and improved by incorporating second order corrections in gradients of the neutron and 
proton densities perturbatively. The calculated symmetry free energy coefficients are found 
to be larger than the symmetry energy coefficients by ~ 10% at medium excitations because 
of the contribution from the symmetry entropy. At low excitations, as expected, there is 
little difference in the values of the two coefficients. At the highest excitation (10 MeV/A) 
that we explore, the difference is ~ 15%. The calculated coefficients Cp compare favorably 
with the available experimental data. 

Both for infinite nuclear matter and for finite nuclei, the symmetry coefficients vary 
linearly with excitation energy; however, for finite systems, the dependence is much stronger. 
The dependence on the EoS of the symmetry coefficients for nuclear matter is found to 
be rather weak. For finite systems the dependence is more noticeable and the symmetry 
coefficients decrease with the softness of the EoS. They are, however, not too sensitive to the 
EoS chosen. The coefficients Ce and Cp of finite nuclei are system dependent, an interplay of 
the role of mass and isospin is quite evident there. For the same asymmetry, the coefficients 
get smaller with smaller mass, this dependence is weak. On the other hand, for the same 
mass, the coefficients show a relatively stronger dependence on the isospin content of the 
nucleus. 

The characterization of the density and excitation dependence of the symmetry term of 
the nuclear interaction is instrumental for the understanding of a plethora of phenomena in 
both nuclear physics and astrophysics. This topic currently attracts much theoretical and 
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experimental activity. Our present calculations have been done in the mean-field framework, 
effects beyond mean-field like clusterization at low densities may have perceptible effects 
and therefore are worth a study. The relevant experimental data are still very scarce, the 
continuing experimental effort in reactions with neutron-rich stable nuclei and future data 
from reactions with exotic isotopes in the radioactive ion beam facilities would contribute 
to a better understanding of these phenomena. 
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Figure Captions 

Fig. 1 The symmetry energy coefficient Ce for nuclear matter at T = as a function of 
density with different interactions. The full line refers to calculations with SBM in- 



teraction; the other results are taken from Ref. 47 



Fig. 2 The symmetry coefficients Ce and Cp for nuclear matter calculated with the SBM 
interaction, (a) shown as a function of density at T = 10 MeV and (6) as a function 
of excitation energy. 

Fig. 3 (Color online) The equilibrium temperature (a), the equilibrium central density (6), 
and the symmetry coefficients (c) as a function of excitation energy for ^^°Sn. The 
calculations are performed with the base density (without self-similar expansion). In 
the bottom panel, the dotted black line and the dash-dot blue line are the symmetry 
coefficients Ce and Cp, respectively, calculated in LDA. The dashed magenta line and 
the solid red line refer to Ce and Cp with inclusion of second order corrections. For 
the experimental data points, see the text. 

Fig. 4 (Color online) Correlation of the symmetry coefficients with density for the system 
^^°Sn. The calculations refer to those with the base density. The experimental points 



are taken from Ref. [ISj]. The notations for the calculated results are the same as in 
Fig. 3(c). 

Fig. 5 (Color online) Same as in Fig. 3, but the calculations are done with the microcanonical 
equilibrium density. 

Fig. 6 (Color online) Same as in Fig. 4, but with the microcanonical equilibrium density. 

Fig. 7 The dependence of the symmetry coefficients on the underlying EoS is shown for the 
finite system ^-^"Sn. 

Fig. 8 (Color online) The symmetry coefficient Ce for the five nuclei studied is shown as a 
function of excitation energy in the upper panel, and as a function of density in the 
lower panel. From top to bottom, the lines correspond to ^^''Sn, ^^^Au, ^^°Sm, ^^°Cs, 
and ^°S, respectively. 
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